o 
o 

(N 

r^ 

o 
< 

d 



> 

00 

m 
o 
l> 

o 

o 
> 

x 

S3 
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The causes of major and rapid transitions observed in biological macroevolution as well as in the 
evolution of social systems are a subject of much debate. Here we identify the proximate causes of 
crashes and recoveries that arise dynamically in a model system in which populations of (molecular) 
species co-evolve with their network of chemical interactions. Crashes are events that involve the 
rapid extinction of many species and recoveries the assimilation of new ones. These are analyzed 
and classified in terms of the structural properties of the network. We find that in the absence of 
large external perturbation, 'innovation' is a major cause of large extinctions and the prime cause of 
recoveries. Another major cause of crashes is the extinction of a 'keystone species'. Different classes 
of causes produce crashes of different characteristic sizes. 



Major transitions in biological and social systems are 
attributed primarily to some combination of external per- 
turbation, selection, novelty, and a complex internal dy- 
namics and structure of the system Jl-pj. While em- 
pirical studies attempt to identify proximate causes for 
individual events (which requires stupendous effort be- 
cause of the complexity of the systems and the difficulty 
of obtaining data) most modeling efforts have tended to 
focus on the statistics of the events taken as a whole 
(see, e.g., [§[|0| for reviews of models in macroevolution 
and finance, respectively). Here we discuss an evolution- 
ary model that permits a precise dissection of the (often 
multiple) causes of individual events. While the model 
is abstract, highly simplified and motivated by chemi- 
cal evolution, the structures and processes that arise in 
it seem to have the flavour of phenomena in biological 
and social evolution. These include the appearance of 
a 'core' and 'periphery' in the network structure of the 
system, a shifting balance between cooperative and com- 
petitive processes as structures evolve, 'core-shifts', his- 
tory dependence, 'keystone species', 'innovations' that 
are 'core-transforming' or 'dormant', and others. Their 
precise mathematical formulation and analysis of their 
role in major transitions in this simplified context may 
help in constructing and analyzing more realistic models. 



The Model 

The system is an idealized prebiotic pond containing a set 
of s chemical species. A given species j can be a catalyst 
for the production of another species i with some small a 
priori probability p. Then the presence of j in the pond 
causes the population of i to increase according to the 
rate equation for catalyzed chemical reactions. The cat- 
alytic relationship is represented graphically by an arrow 
from node j to node i in a directed graph representing 
the chemical network. The graph is completely specified 
by its adjacency matrix C — (cy), i,j € {1, . . . , s}. cy 
is unity if there is a link from node j to node i, i.e., if 
species j catalyzes the production of species i, and zero 



otherwise. Each species has a population y,i in the pond. 
Xi = Hi/ J2j Vj i s ^s relative population; by definition 

< Xi < 1, X)i=i Xi = I. On a certain timescale T the 
relative populations of the species reach a steady state 
denoted X = (Xi, . . . , X s ) that depends on the catalytic 
network. We imagine the pond to be subject to peri- 
odic external perturbations on a timescale greater than 
T in the form of tides, storms or floods. Such a pertur- 
bation can flush out existing molecular species from the 
pond (the ones with the least Xi being the most likely 
to be eliminated) and bring in new chemical species from 
the environment whose catalytic relationships with ex- 
isting species in the pond are quite different from the 
eliminated ones. Then, after the perturbation, the popu- 
lations in the pond will evolve to a new steady state, be 
subject to another perturbation, and so on. 

The precise rules are as follows: Initialization: Start 
with a random graph with s nodes and 'catalytic proba- 
bility' p. That is, for every ordered pair (i,j) with i =/= j, 
Cij is unity with probability p and zero with probability 

1 — p. cu = for all i (to forbid self- replicating species). 
Fig. la is an example. Assign each x\ a random num- 
ber between and 1 and uniformly rescale all Xi so that 
Si=i x i ~ 1- Dynamics: First, keeping C fixed, evolve 
x from its initial condition according to 
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(1) 



for a time T large enough to reach its attractor. De- 
note Xi (T) = Xi . Eq. [j] follows from the rate equation 
Vi = Y2j CijVj ~ 4>Di f° r th° populations, which in turn is 
an idealization of rate equations in a well stirred chemical 



reactor. 1 Find the set of nodes with the least Xi. Sec- 
ond, pick a node (denoted k) from this set at random and 
remove this node from the graph along with all its links. 
Add a new node (also denoted k) to the graph which is 
connected randomly to the existing nodes according to 
the same catalytic probability p. Mathematically, this 
means that for every i ^ fc, Cki and Cik are reassigned to 
unity with probability p and zero with probability 1 — p, 
irrespective of the value they had earlier, and c^u = 0. 
Set Xk = Xq (a small constant) , perturb all other Xi about 
their existing value X t by a small amount, and uniformly 
rescale all Xi to preserve the normalization X)i=i Xi = 1. 
This procedure provides a new graph and a new initial 
condition for x. Now return to the first step of the dy- 
namics and iterate the procedure several times. 2 

This model, introduced in rcf. p!l| , was inspired by the 
work in refs. jj^-[L5Pf| . The removal of the least pop- 
ulated species implements selection J5J, and its replace- 
ment by another randomly connected species implements 
the introduction of novelty into the system. The network 
character is shaped by the repeated application of these 
two incremental external perturbations along with the 
internal population dynamics. 

Three regimes of behaviour and transitions be- 
tween them. The system exhibits three regimes or 
phases of behaviour. This is illustrated in Fig. 2 which 
shows the number of populated species in the attractor, 



S\i (those species for which Xi > 0) vs. time (n, the 
number of graph updates) for a run with s = 100 and 
p = 0.0025. In the 'random phase' s± stays low with small 
fluctuations. In the 'growth phase', s± typically rises ex- 
ponentially with occasional drops. Finally, in the 'orga- 
nized phase' si stays close to s, the maximum value it 
can take. The random and growth phases were discussed 
in (^J^,[17|]. As is evident from Fig. 2, the organized 
and the growth phases exhibit occasional discontinuous 
transitions or 'crashes' in which a number of species sud- 
denly go extinct (their Xi become zero in a single time 
step). At the end of a crash the system is in the random 
or growth phase. This is followed by a recovery in which 
the system moves again towards the organized phase. In 
|18| it is shown that crashes are primarily 'core-shifts', 
a specific kind of change in the structure of the graphs 
(discussed below) and recoveries are due to the growth 
of 'autocatalytic sets'. The main purpose of this paper 
is to elucidate the mechanisms which cause core-shifts. 
Crashes and recoveries also occur in a related model with 
negative links and variable link strengths fll7| ] , where they 
are more difficult to study analytically. 

Definitions and notation 

Autocatalytic set (ACS): An ACS is a subgraph each of 
whose nodes has at least one incoming link from a node 
of the same subgraph. (By a subgraph we mean a subset 
of nodes together with all their mutual links.) Thus an 



x The rate equation ijt = fc(l + vyj)nAns — 4>yi follows from 
the reaction scheme A + B —* i, where A and B are reactants 
with populations ua and t%b, j and i are catalyst and prod- 
uct with populations t/j and yi respectively, and is a death 
rate or dilution flux in the reactor (k is the rate constant for 
the spontaneous reaction and v is the catalytic efficiency). 
If riA,nB are large and fixed and the spontaneous reaction 
is much slower then the catalyzed reaction, this equation re- 
duces to yi = cyj — cpyi, where c is a constant. A generalization 
of the latter equation is yi — X^=i Ci iVi ~ $]& f° r the case 
where species i has multiple catalysts. Eq. Fy follows from 
this by taking the time derivative of Xi — yi/^2 ._, 2/j. In 
the present model we make the idealization that all catalytic 
strengths are equal. The second (quadratic) term in Eq. Fy 
is needed to preserve the normalization of the Xi under time 
evolution. Note that it follows automatically from the nonlin- 
ear relationship between Xi and yi when the time derivative 
of Xi is taken. 

2 The attractor configuration X is determined in this paper 
by its algebraic properties discussed later, not by numerically 
integrating Eq. hi Hence we are effectively taking T — oo 
which allows us to prove certain results analytically. Prelim- 
inary simulations suggest that when Eq. hi is numerically in- 
tegrated for a fixed, but large, T, qualitatively similar results 
are obtained, although one can expect that some of the sharp 
transitions are now spread out over several graph update time 
steps. 



ACS contains a catalyst for each of its members [[19] £1| . 
In each of Figs, lb-k, the subgraphs formed by the set of 
all black nodes or all the black and grey nodes are ACSs. 
Figs, la and 11 do not have an ACS. For any subgraph A, 
let Xi(A) be the largest eigenvalue of the adjacency ma- 
trix of A. We denote Ai(C) = Ai. It can be shown ]lq] 
that if the graph does not have an ACS, then Ai = 0, 
and if it does, then Ai > 1. Ai therefore represents a 
topological property of the network. 
Dominant ACS: It can be shown that if Ai > 1 the sub- 
graph comprising the populated species (X t > 0) must 
be an ACS |ll],|l6| , which will be referred to as the domi- 
nant ACS. The dominant ACS is uniquely determined by 
the graph and does not depend upon the initial condition 
for x (except for special initial conditions forming a set 
of measure zero, which we ignore). The subgraph formed 
by the set of all black and grey nodes in each of Figs. 
lb-k is the dominant ACS for that graph. In addition to 
its topological significance, Ai also has a dynamical in- 
terpretation as being the 'population growth rate' of the 
dominant ACS. 3 

Core and periphery of the dominant ACS: The core of 
the dominant ACS of a graph C (sometimes also referred 
to as the 'core of C") is the maximal subgraph, Q, from 
each of whose nodes all nodes of the dominant ACS can 
be reached along some directed path. The rest of the 
dominant ACS is its periphery. The subgraph of all black 
nodes in Figs, lb-k constitutes the core of the dominant 
ACS in the graph. 4 Every periphery node has an incom- 
ing path that originates from the core, but no outgoing 
path that leads to the core. The periphery can contain 
loops within itself (e.g., the 2-cycle between nodes 36 and 
74 in Fig. le). One can prove that Ai(Q) = Ai(C) = Ai. 
A subgraph with more than one node containing a di- 
rected path from each of its nodes to each of its other 
nodes is said to be irreducible 22 1. An irreducible sub- 
graph is always an ACS but the converse is not true. 
However, the core (of each component of a dominant 
ACS) is an irreducible subgraph. Since Xi(A) > 1 for 
any irreducible subgraph A, it follows that the latter is 
a 'self-sustaining' structure in the sense that if no other 
links were present in the graph, the nodes of A would 



still have nonzero AQ by virtue of their mutual links. Ai 
measures the 'strength' of the core in two ways: one, its 
intrinsic population growth rate, and two, its multiplicity 
of internal pathways. To see the latter, compare the in- 
creasing and decreasing pattern of Ai between n = 2854 
and 5042 in the inset of Fig. 2 with the sequence of Figs, 
lb-f. When the core (of every disjoint component of the 
dominant ACS) has exactly one cycle (Figs. lb,f-i,k), 
then Ai = 1, and vice versa. Such a core is fragile due to 
the absence of any redundancy in its internal pathways; 
the removal of any link from such a core will cause the 
ACS property (of that component) to disappear. 
Crash: a crash is a graph update event, n, in which a 
significant number (say > s/2) of the species go extinct, 
i.e., Asi(n) = Si(n) — si(n — 1) < —s/2. 
Core overlap: given any two graphs C and C whose 
nodes are labeled, the core overlap between them, de- 
noted Ov(C,C), is the number of common links in the 
cores of C and C . If any one of them does not have an 
ACS, Ov(C,C') = 0. 

Core-shift: A core-shift is a graph update event in which 
Ov(C n -i,C n ) = 0, i.e, there is no overlap between the 
cores of the dominant ACS before and after the event 
(C„ is the graph at time step n). 

Keystone species. One can consider the impact of the 
hypothetical removal of any species i from the graph, ir- 
respective of whether i is the least populated or not. For 
example one can ask for the core of the graph C — i that 
would result if species i (along with all its links) were 
removed from C . A species i will be referred to as a 
keystone species if Ov(C, C — i) = 0. Thus a keystone 
species is one whose removal modifies the organizational 
structure of the graph (as represented by its core) dras- 
tically. In the ecological literature a keystone species 
is regarded as one whose elimination from the ecosys- 
tem would cause a significant fraction of species in the 
ecosystem to go extinct p3H26fl. We will see that like- 
wise in the present model the removal of keystone species 
causes large crashes. (Indices of keystoneness based on 
Ov(C,C — i) or on the change in s\ caused by the re- 
moval of a species can also be defined.) Note that if 



This follows from the fact that the attractor configura- 
tion X is always an eigenvector of C with eigenvalue Ai, 
i.e., ^.CijXj — \\Xi fill]. Thus, when (j> — 0, substitut- 
ing yi (X Xi in the population dynamics equation y = Cy, 
one gets y = Aiy. 

4 Sometimes the dominant ACS consists of two or more dis- 
joint subgraphs, as in Fig. li. Then the above definition 
applies to each component separately. There exist other ACS 
structures for which this definition is not adequate, e.g., two 
disjoint 2-cycles pointing to a single downstream node. Such 
structures arise rarely and can be treated by a more general 
definition of core and periphery without altering the main 
conclusions presented here. 



Ai = 1 and the ACS has a single connected component 
(e.g., figs. lb,f-h,k), every core node is necessarily a key- 
stone species since its removal would destroy the cycle 
that constitutes the core. 



Innovations. In the present model the new species k 
at each time step together with its set of new links may 
be regarded as a 'novelty' introduced into the system. In 
the new attractor the new species may go extinct, i.e., X^ 
may be zero, or it may survive, i.e., X^ is non-zero. Let us 
define an innovation to be a novelty in which the relative 
population of the new species in the new attractor just af- 
ter the novelty occurs is nonzero. This definition has the 
feature that an innovation always involves new connec- 
tions. It does not use any exogenously defined notion of 
fitness. The only performance criterion it requires is that 
the new links should enable the new node to survive until 
the next graph update. Even this minimal requirement 
has nontrivial consequences. For instance when the new 
species receives an incoming link from an existing dom- 
inant ACS, it typically has a nonzero population in the 
new attractor. Each recovery process, which occurs due 
to the expansion of the dominant ACS during the growth 
phase 1 18 1, is an accumulation of just such innovations. 

We will be interested in a special class of innovations 
in which the novelty creates a new populated irreducible 
subgraph. Such innovations create or add to a 'self- 
sustaining structure' in the graph, in the sense discussed 
earlier. For instance the appearance of the first ACS at 
n = 2854 (see Fig. lb) is such an innovation. There, 
species 90, which was a singleton before the event, went 
extinct and was replaced by a new species 90 that had an 
incoming and outgoing link to node 26. The two formed 
an irreducible subgraph whose Ai value was 1. This in- 
novation in fact triggered the self organization of the net- 
work around this ACS. By n — 3880 (see Fig. lc) the 
core had grown to 18 nodes as a result of several events in 
which the new species was an addition to the core. Every 
event in which the core is strengthened by the addition 
of a new node is also an innovation in which a new pop- 
ulated irreducible subgraph is formed. When such an 
innovation arises, denote the new irreducible subgraph 
which includes the new species as N (or N n at time step 
n.) N or N n will stand for the maximal irreducible sub- 
graph of which the new species k is a member. It fol- 
lows from the Perron-Frobenius theorem [E2J that such 
an innovation necessarily increases the Ai value of some 



substructure in the graph, i.e., \i(N) > Ai(iV — k). N n 
becomes the new core of the graph if it is 'stronger' than 
the old core. More precisely, one can show that N n will 
become the new core of the graph, replacing the old core 
Q n -i, whenever any one of the following two conditions 
hold: 

(a) Ai(JV n ) > XxiQ'J or, 

(b) Ai(A„) = Xi(Q' n ) and N n is 'downstream' of Q' n . 5 
Such an innovation may be referred to as a core- 
transforming innovation. If Q n -\ C N n , such an inno- 
vation enlarges the existing core. However, if Q n -\ and 
N n are disjoint, we get a core-shift. 

In the present model, the generation of novelty does 
not depend upon existing structure since the links of a 
new node are chosen randomly from a fixed probability 
distribution: novelty is 'noise' in this model. (Variants of 
the model that depart from this are easily constructed.) 
However, whether the novelty constitutes an innovation 
is 'context dependent' (i.e., dependent upon the struc- 
ture of the existing network). Also, the short and long 
term impact of an innovation depends upon the (histor- 
ical evolution of the) 'context', as will be seen below. 



Classification of core-shifts 

In a set of runs with s — 100, p — 0.0025 totaling 1.55 
million iterations we observed 701 crashes. 612 of these 
were core-shifts flq]. Fig. 3 differentiates between the 
612 core-shifts we observed. They fall into three cate- 
gories: (i) complete crashes (136 events), (ii) takeovers 
by core-transforming innovations (241 events), and (iii) 
takeovers by dormant innovations (235 events). 



Complete crashes. A complete crash is an event in 
which an ACS exists before but not after the graph up- 
date. Therefore for a complete crash at time step n, 
Ai(C n _i) > 0andAi(C„) = 0. These events take the sys- 
tem back to the random phase. For example at n = 8232, 
in Fig. Ik, node 54 is one of the least populated species 
and is hit at n = 8233. It is replaced by a new species 
that has a single outgoing link to node 50 and no in- 
coming link, resulting in the complete disruption of the 
ACS. It is evident that complete crashes must always be 
caused by the elimination of a keystone species. Further- 
more, Fig. 3 shows that Ai(C n _i) = 1 for every complete 
crash observed in the runs. Hence the core of the ACS 
is a single cycle when the event occurs and the species 



5 We use the notation C' n = C n -i — k for the graph of s — 1 
nodes just before the novelty at time step n is brought in 
(and just after the least populated species k is removed from 
C„_i). Q' n stands for the core of C' n . A subgraph A is 'down- 
stream' of another subgraph B if there exists a directed path 
from some node of B to a node of A but none from any node 
of A to a node of B. 



removed is a member of that cycle. 

Takeovers by core-transforming innovations. An 

example of a takeover by a core-transforming innovation 
is given in Figs. lg,h. At n — 6061 the core was a 
single loop comprising nodes 36 and 74. Node 60 was 
replaced by a new species at n = 6062. The new node 
60 created an innovation at n = 6062, with Neo62 being 
the cycle comprising nodes 60,21,41,19 and 73, down- 
stream from the old core. The graph at n = 6062 has 
one cycle feeding into a second cycle that is downstream 
from it. (This is an example of condition (b) for a core- 
transforming innovation.) The population attractor of 
such a graph has the property that the Xj of the nodes 
in the upstream cycle are all zero and only the nodes of 
the second cycle (as well as nodes further downstream 
from it) have nonzero Xi. Thus when the above inno- 
vation arises, the new cycle becomes the new core and 
all nodes that are not downstream from it get Xi — 0, 
resulting in a large drop in s\ from 89 to 32. For all such 
events in Fig. 3, \i(Q' n ) = Ai(C n _i) since k happened 
not to be a core node of C n -i- Thus these core-shifts 
satisfy Ai(Cn) = \i(N n ) > Xi(Q' n ) = Ai(C„_i) > 1 in 
Fig. 3. 

Takeovers by dormant innovations. Figs. le,f show 
an example of a takeover by a dormant innovation. At 
n = 5041, the core has Ai = 1.24 and there is a cycle com- 
prising nodes 36 and 74 in its periphery. Node 85 is hit 
which results in a cycle (26 and 90) feeding into another 
cycle (36 and 74) at n = 5042. Thus at n = 5042, for the 
same reason as in the previous paragraph, 36 and 74 form 
the new core with only one other downstream node, 11, 
being populated. All other nodes become depopulated 
resulting in a drop in s\ by 97. Such a core-shift is the 
result of an innovation that arose earlier (the cycle be- 
tween 36 and 74 arose at n = 4695) but lay dormant 
downstream of the existing core until one of the keystone 
species of the latter (node 85) was hit and made it weak 
(i.e., reduced its Ai to a value less than or equal to the 
Ai of the downstream innovation). A dormant innovation 
can takeover as the new core only following a keystone 
extinction which weakens the old core. In such an event 
the new core necessarily has a lower (but nonzero) Ai 
than the old core, i.e., Ai(C„-i) > \\{C n ) > 1. Note 
that at n = 5041 if the downstream cycle between 36 
and 74 were absent, 85 would not be a keystone species 
by our definition, since its removal would still leave part 
of the core intact (nodes 26 and 90). 85 becomes key- 
stone, and the core of which it is a part becomes fragile 
and susceptible to a core-shift, because a self-sustaining 



innovation has occurred in the distant periphery. 

If we consider the set of all drops in. si, large or small, 
where an ACS exists before the event (there are 126454 
such events in the above mentioned runs) we find that the 
number of complete crashes remains the same, 136 (the 
mean size of the drop in such events is | Asi | = 98.2 with 
a standard deviation a = 1.2), while core-shifts caused 
by dormant innovations go up to 359 (with |Asi = 62.2, 
a = 25.9) and those due to core-transforming innova- 
tions to 524 (|Asi| = 48.2, a = 25.6). The rest of 
the events consist of 9851 'partial core-shifts' (in which 
the core changes, but Ov{C n -i,C n ) ^ 0; in this cate- 
gory |Asi| = 2.18, a = 7.42), and 115584 events where 
there was no change in the core but the periphery is af- 
fected (|Asi| = 1.05, a = 0.99). Thus different classes of 
proximate causes of drops arise dynamically with differ- 
ent frequencies and typically produce events in different 
size ranges. The ranges, however, overlap and some dis- 
tributions have fat tails (e.g., of the 701 crashes with 
|Asi| > 50, there are 79 and 10 events, respectively, in 
the last two categories) . Detailed distributions and their 
dependence on s,p are open questions (but see jig] )- 

Discussion 

The present model exhibits mechanisms by which inno- 
vations can play a major role in crashes and recoveries 
in a complex system. 6 It provides a mathematical ex- 
ample of 'creative destruction' |27| at work in causing 
large upheavals. It distinguishes two processes involv- 
ing innovations, both having analogues in the real world. 
One is exemplified by the appearance of the automobile 
which made the horse drawn carriage and its ancillary 
industries obsolete. This is like the example of the core- 
transforming innovation shown in Figs. lg,h where the 
graph update produced a self-sustaining structure that 
was more vibrant than the existing core within the con- 
text of the present organization. This structure became 
the new core, rendering many nodes drawing sustenance 
from the old core dysfunctional. The subsequent devel- 
opment of other industries dependent on the automobile 
mirrors the growth of the ACS around the new core. The 
second process is exemplified by the emergence of the 
body plans of several phyla which are dominant today. 
It is believed that while these body plans originated in 
the Cambrian era more than 520 million years ago P9[ , 
the organisms with these body plans played no major role 
till about 250 million years ago. They started flourish- 
ing only when the Permian extinction depleted the other 
species that were dominant till that time JMJ. This is 
similar to the events shown in Figs. le,f where an earlier 
innovation had lain dormant for a while without disturb- 



6 Analogues of innovations and core-shifts appear to be play- 
ing an important role in another related but quite different 
model Pq] where rapid transitions are observed. 



ing the existing core, but when the latter became suffi- 
ciently weak, took over as the new core and flourished. 

Recently there has been substantial progress in graph 
theoretic analyses of complex systems and in particu- 
lar 'small-world' |3l] and 'scale-free' |32| properties have 
been found for several real networks (for reviews and 
references see |33|,|34|]). It may be interesting to study 
whether certain classes of real networks also have some 
kind of a 'core-periphery' structure, in view of its possible 
dynamical significance. 
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Figure legends 

Figure 1. The structure of the evolving graph at various 
time instants for a run with s = 100, p = 0.0025. Node 
numbers i from 1 to 100 are shown in the circles repre- 
senting the nodes. Nodes with zero relative population 
in the attractor configuration for the graph (Xi = 0) 
are shown as white circles; the rest (si in number) have 
nonzero Xi. In the graphs where an autocatalytic set 
(ACS) exists, black circles correspond to nodes in the 
'core' of the ACS, and grey to the 'periphery', defined 
in the text, (a) n = 1, the initial random graph, (b) 
n = 2854, where the first ACS, a 2-cycle between nodes 
26 and 90, appeared, (c) n = 3880, the beginning of the 
organized phase when the ACS first spanned the entire 
graph, (d) n = 4448, when the core reached a peak in 
the number of loops it contained, (e) n — 5041, just be- 
fore a 'core-shift', (f) n = 5042, just after the core-shift 
caused by a 'keystone' extinction in the presence of a 
'dormant innovation', (g) n — 6061, just before another 
core-shift, (h) n = 6062, just after the core-shift caused 
by a 'core-transforming innovation', (i) n — 6070, when 
the old core stages a come back as a disconnected com- 
ponent after node 32 becomes a singleton, (j) n = 6212, 
when the new core strengthens itself and depopulates the 
recently resurgent old core, (k) n = 8232, just before the 



first 'complete crash', (1) n — 10000, between the first 
'complete crash' and the subsequent 'recovery'. 

Figure 2. The number of populated species, s\ (contin- 
uous line), and the largest eigenvalue of C (whose sig- 
nificance is discussed later in the text), Ai (dotted line), 
versus time, n for a run with s = 100 and p = 0.0025. 
The Ai values shown are 100 times the actual Ai value. 
The first 10000 time steps are enlarged in the inset. The 
run is the same for which the graph snapshots arc shown 
in Fig. 1. The impact of the events described in Fig. 1 is 
clearly visible in this curve. E.g., at n — 2854 Ai jumps 
from zero to one and si exhibits the first sustained up- 
ward movement, at n — 3880 Si hits its maximum value, 
100, and then fluctuates mainly between 99 and 100, and 
at n = 4448 Ai reaches a local maximum, si drops from 
100 to 3 as a result of the 'core shift' at n = 5042, and 



from 89 to 32 in the core-shift at 6062. At n = 6070 a 
large recovery event is seen as the old core and the still 
intact part of its periphery get repopulated, only to be 
trounced again at n — 6212 when the new core strength- 
ens itself to a Ai value greater than 1. At n = 8233 
Si crashes from 100 to 2 when the ACS is completely 
destroyed and Ai drops from one to zero. 

Figure 3. Frequency, /, of the 612 core-shifts ob- 
served in a set of runs with s = 100 and p = 0.0025 
vs. the Ai values before, Ai(C„_i), and after, Ai(C n ), 
the core-shift. Complete crashes (black; Ai(C„_i) = 1, 
Ai(C n ) = 0), takeovers by core-transforming innovations 
(blue; Ai(C„) > Ai(C„-i) > 1) and takeovers by dor- 
mant innovations (red; Ai(C„-i) > Ai(C n ) > 1) are dis- 
tinguished. Numbers alongside vertical lines represent 
the corresponding / value. 
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